Pairing Correlations and the Pseudo-Gap State: Application of the "Pairing 

Approximation" Theory 
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We investigate the pseudogap onset temperature T* , the superconducting transition temperature Tc 
and the general nature of the pseudogap phase using a diagrammatic BCS-Bose Einstein crossover 
theory. This decoupling scheme is based on the "pairing approximation" of Kadanoff and Martin, 
further extended by Patton (KMP). Our consideration of the KMP "pairing approximation" is driven 
by the objective to obtain BCS like behavior at weak coupling, (which does not necessarily follow for 
other diagrammatic schemes). Two coupled equations, corresponding to those for the single particle 
and pair propagators, must be solved numerically, along with the number equation constraint. The 
variation of small to large coupling constant g is explored, whereby the system is found to cross 
over from BCS to Bose-Einstein behavior. Our numerical calculations proceed in two stages: first 
we investigate the "lowest order theory", which is appropriate at temperatures well above Tc- We 
use this theory to determine where the Fermi liquid state first breaks down. This breakdown, which 
occurs at T* and is associated with intermediate values of the coupling, corresponds to a splitting of 
the single peaked (Fermi liquid) electronic spectral function into two peaks well separated by a gap, 
as might be expected for the pseudogap phase. Indeed, our calculations provide physical insight into 
the pseudogap state which is characterized by the presence of metastable pairs or "resonances" , which 
occupy states around the Fermi energy; in this way they effectively reduce the single particle density 
of states. The superconducting instability Tc is evaluated in the second stage of our calculations. 
Here we introduce "mode coupling" efi'ects, in which the long lived pairs are affected by the single 
particle pseudogap states and vice versa. Our Tc equations, which turn out to be rather simple as 
a result of the KMP scheme, reveal a rich structure as a function of g in which the pseudogap is 
found to compete with superconductivity. Our results are compared with alternate theories in the 
literature. 



PACS numbers: 74.20.Mn, 74.25.-q,74.25.Fy, 74.25.Nf, 74.72.-h 



I. INTRODUCTION 

Pseudogap phenomena were reported in the first gener- 
ation of underdoped high Tc samples [Q, although they 
received little attention at that time. These early ex- 
periments 01 on underdoped cuprates demonstrated a 
suppression in the magnetic susceptibihty with decreas- 
ing temperature (T < T*) which was suggestive of some 
form of singlet pairing in the normal state. These first in- 
dications of the pseudogap state were strictly associated 
with magnetic probes. Hence the phrase "spin gap" was 
used to describe them and related frequency dependent 
magnetic phenomena ||]. More recently, it has become 
clear that similar effects Q are present in the charge 
response as well. Today the spin gap is more broadly 
interpreted to correspond to some form of excitation gap 
(i.e., pseudogap) in both charge and spin channels. It 
can and has been debated as to whether the effects are 
equally profound in each. 

Among the most revealing data on the pseudogap 
state are angle resolved photoemission experiments [Q]^] 
(ARPES) which demonstrate that below T* the pseu- 
dogap is present in the energy dispersion of the normal 
state electrons and that, moreover, it seems to have the 
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same (d-wave) symmetry as the superconducting state, 
and evolves more or less continuously through Tc into a 
superconducting excitation gap. Finally, these data show 
that the pseudogap is attached to an underlying Lut- 
tinger volume Fermi surface, although there are many 
aspects about the high temperature regime (T > T*), 
(where the pseudogap is essentially absent), which do 
not correspond to a canonical Fermi liquid. 

Theories of the cuprate pseudogap can be broadly di- 
vided into two groups: those that associate this state 
with magnetic pairing of some sort [^-0 and those which 
view the pseudogap as deriving from some form of precur- 
sor superconductivity [p^-p2 24 1. (In reality the division 
between these two classes is somewhat blurred, since in 
many scenarios magnetic pairing is a natural requisite for 
superconductivity). The present paper belongs to this 
second category. It is our contention that this second 
scenario deserves careful attention in large part because 
mean field descriptions of the transition (which neglect 
precursor or pairing correlation effects) are not expected 
to be applicable to short coherence length, quasi 2d su- 
perconductors such as the cuprates. Moreover, ARPES 
data suggest a smooth evolution of the pseudo- into the 
superconducting gap. Finally, these systems ultimately 
undergo a superconducting, not magnetic transition. At 
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the very least it is clear that the pseudogap state will 
never be fully characterized unless precursor supercon- 
ducting effects are properly calibrated. 

Within the precursor superconductivity school, there 
are, moreover, two distinct classes. The first of these is 
based on the short coherence length ^ which is naturally 
associated with modest size Cooper pairs and, therefore, 
a superconducting state intermediate between the BCS 
(large pair) and Bose-Einstein (small pair) descriptions. 
The second class is based on the anomalously low plasma 
frequency ujp which suggests that phase fluctuations play 
an important role. In this way the amplitude of the order 
parameter is established at T*, while phase coherence 
occurs at the much lower temperature Tc- In reality, 
both ^ and ujp are small and some initial work by our 
group [p5| has treated these two parameters on a more 
even footing. 

While, there have been no decisive arguments in favor 
of one or the other of the two precursor superconductivity 
schools, in this paper we explore the "BCS Bose-Einstein 
crossover" scenario, in large part because it lends itself 
to more detailed diagrammatic and therefore quantita- 
tive calculations and insights. There is a long history 
associated with this small pair theory. Leggett first 
proposed a variational ground state wave function which 
exhibited a smooth interpolation between the BCS and 
Bose-Einstein regimes. An essential component of his 
approach was the introduction of a self consistent equa- 
tion for the chemical potential to be solved in conjunc- 
tion with the variational conditions [ p7[ . This fi differs 
significantly from the Fermi energy in all but the very 
weak coupling limit. Nozieres and Schmitt-Rink (NSR) 
p8[ extended Leggett's ideas to non-zero temperature 
and thereby deduced a transition temperature Tc which 
varied continuously with increasing coupling constant g 
from the BCS exponential dependence on g to the Bose- 
Einstein asymptote. Randeria and co- workers fl^ were 
among the first to apply the NSR approach, derived via a 
saddle point scheme, to the cuprates, which were associ- 
ated with intermediate values of g. Subsequent extensive 
numerical ||2^ and analytical studies have addressed 
the intermediate g regime beyond the NSR approxima- 
tions and find a variety of aspects of pseudogap behavior. 

An important goal of the present paper is to address 
the cross-over problem within a specific diagrammatic 
scheme, which is amenable to quasi-analytic calculations 
and physical interpretation. In the process we are able to 
compute the electronic self energy, spectral function and 
a quantitative phase diagram for Tc and T* as a function 
of g. We refer to this decoupling scheme as the "pair- 
ing approximation" theory. It is based on early work 
by Kadanoff and Martin and has been argued by 
these authors to be the most appropriate treatment for 
introducing pairing correlations of the kind encountered 
in a conventional superconducting state. Moreover, this 
approximation has been extensively applied by Patton 



p2| , ^ to the problem of T > Tc fluctuations in low di- 
mensional, dirty superconductors. As will be discussed 
in detail in a future paper |Q| which addresses the su- 
perconducting state, the present diagram scheme leads 
naturally to Leggett's |^ description of the interpolated 
ground state. This is in contrast to alternative diagram- 
matic approaches ]35|-|39|] which, below Tc, do not recover 
the BCS limit in the weak coupling regime. 

As a result of these calculations, a simple physical pic- 
ture of the pseudogap phase emerges. This is a regime of 
intermediate coupling g which is characterized by a nor- 
mal state that is somewhere between the "free" fermions 
of the BCS limit and the "pre-formed" pairs of the Bose- 
Einstein case. A natural description of the interpolated 
state is to presume that it corresponds to "resonantly" 
scattered or meta-stable pairs. Indeed, our calculations 
of the electronic self energy and spectral function ]2^ ] 
show that the presence of resonant scattering (as deduced 
from the character of the pair propagator or T-matrix), 
correlates rather well with the break-down of the Fermi 
liquid phase. Once the pairs are meta- stable they block 
available states around the Fermi surface which would 
otherwise be occupied by single electrons and thereby 
introduce a pseudogap into the electronic spectrum. 

While this interpretation of the density of states de- 
pression is in the same spirit as earlier work on super- 
conducting fluctuation induced spectral gaps, the present 
description differs in an important way. Long lived pair 
states in ordinary superconductors occur in the vicinity 
of Tc as a consequence of critical slowing down. By con- 
trast, here they are found to occur at much higher tem- 
peratures simply because the coupling g is sufficiently 
strong to bind the pairs into a meta-stable state. This 
picture must necessarily be distinguished from the phase 
fluctuation picture of Emery and Kivelson jl2) as well 
as from the "pre-formed" pair models of other groups 
|l5|-p^. Our pairs correspond to microscopic as distinct 
from mesoscopically established regions of superconduc- 
tivity. Moreover, the resonant pairs of the present model 
have a finite lifetime and spatial extent and do not obey 
Bose statistics. 

II. THEORETICAL FRAMEWORK 

We consider a generic Hamiltonian consisting of 
fermions in a three dimensional, jellium gas in the pres- 
ence of an s-wave attractive interaction, Vk,k' — 5"^k¥'k', 
where ipj^ — {1 + fc^/fcg)~^/^ and g < is the coupling 
strength expressed in units of gc = — 47r/mfco- Effects as- 
sociated with quasi two dimensional lattices and d-wave 
pairing have also been addressed within this formalism 
and will be discussed in a companion paper pO[ . The 
results which we present here are generally robust, al- 
though the relevant energy scales decrease considerably 
as the system becomes more two dimensional. Because 
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we consider only this simplest case, we defer the discus- 
sion of direct comparison with ARPES data to future 
work. 

Following earlier work |3|,Q by Kadanoff and Martin 
and by Patton, the system can be characterized by the 
integral equations of the "pairing approximation" which 
describe a coupling between the electronic self energy S 
and T matrix t. We add to these equations the usual 
expression for particle number conservation 
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where the Green's function is given by G 

^Ktuj,^ the bare propagator is G^^^'^ 
^m/^i are the even/odd Matsubara frequencies, with 
electronic dispersion given by ek = k'^ /2m — /i. 

These first two equations may be readily written dia- 
grammatically as shown in Figure |l|. For completeness 
we also indicate the diagrams of the "lowest order the- 
ory". Figure |l|(b), which will be used initially to build 
intuition. 

Although, these integral equations must in general be 
solved numerically, it will be possible to arrive at consid- 
erable analytic insight. The goal of this paper is to solve 
these three coupled equations as a function of arbitrary 
coupling g for two different temperature regimes. Our 
calculations will accordingly proceed in two stages. First 
we discuss the "lowest order" conserving theory (Section 
HI) where all Green's functions in Eqs.(|l|) are replaced 
by bare propagators j2^] . This theory is adequate for de- 
scribing the break down of the Fermi liquid at T* and for 
higher temperatures. In the next stage we (Section |^) 
introduce feedback or "mode coupling" effects in which 
the pair propagator is allowed to depend on the single 
particle self energy and vice versa. This feedback is nec- 
essarily important within the pseudogap phase and all 
the way down to Tc. For a range of temperatures at and 
above Tc but less than T* , where the pseudogap is still 
well established, feedback or mode coupling effects can 
be readily introduced via a simple parameterization of 
the self energy. We will see below that this simplicity oc- 
curs in large part because of the presence of one "bare" 
Green's function in the first two equations above. This 
parameterization will be justified numerically in consid- 
erable detail. By contrast, a direct numerical attack on 
the problem using the three full equations is difficult to 
implement as well as assess because of restrictions to a 
finite number of Matsubara frequencies ||2^. Moreover, 
it does not provide the same degree of physical insight, 
nor can it be done in as controlled a fashion. 



The present scheme should be contrasted with two al- 
ternative approaches in the literature. The approach of 
Nozieres and Schmitt- Rink Q is embedded ||3^] in the 
limit of Eqs. (^ which corresponds to the lowest order 
theory and in which also the right hand side of the num- 
ber equation is approximated by the first two terms in a 
Dyson expansion 
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As a consequence of this last approximation there are 
general problems with the violation of conservation laws 
which can be immediately corrected by summing the full 
Dyson equation for G. In this way our conserving "low- 
est" order theory is obtained as a first stage improvement 
over the NSR scheme. Along related lines, the Baym- 
Kadanoff criteria have been applied |4^,^^ to verify con- 
servation laws in the context of the full, mode coupled 
theory in Eqs. Finally, Ward identities can be im- 
posed in conjunction with these full equations to arrive 
at a satisfactory description of two particle properties as 
has been extensively investigated in Ref. [ p2[ . 

Much attention has been paid to an alternative scheme 
in which all bare Greens functions in Eqs. (|^) are re- 
placed by their dressed counterparts. This approach, 
often referred to as the fluctuation exchange scheme or 
FLEX, has been studied by direct numerical methods in 
finite size systems 39 1 as well as for the purposes of com- 
puting Tc [^. While it is a derivable" theory in the 



sense of Baym 1 42 it was noted some time ago |3l| to be 
inconsistent with BGS theory, even in the weak coupling 
limit. It is, consequently, useful to explore the present 
Kadanoff-Baym-Patton approach both as a point of com- 
parison for results obtained within the FLEX scheme, 
and because it appears to be more analytically tractable 
and more readily interpreted within the context of con- 
ventional theories of superconductivity. 



III. THE LOWEST CONSERVING ORDER 

The single particle self-energy, Sk.w , at the lowest con- 
serving level of approximation corresponds to 
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q,n„ 

and the T-matrix is given by 
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These equations are indicated diagrammatically in Fig- 
ure 0(b). 

Although this formulation of the problem will not give 
quantitatively correct values for properties such as the 
transition temperature, it does yield qualitatively cor- 
rect results for the self-energies and spectral functions. 
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Indeed pseudogap effects in the spectral properties of the 
system are produced already at this order |^2|. This ap- 
proach becomes increasingly more reliable in the weak 
pseudogap regime, where feedback effects are expected 
to be neghgible. 

In the next two subsections we discuss and plot the 
behavior of the T-matrix t^^\^ and, with these insights, 

then address the electronic self energy S^"^ which is de- 
termined from an integral over the T-matrix. Fundamen- 
tal to our considerations is the observation that for suffi- 
ciently strong coupling, when the fermions are bound into 
composite bosons, a Fermi liquid description is clearly in- 
appropriate. The goal of the present discussion is to de- 
termine for what values of g does this breakdown of the 
(normal state) Fermi liquid happen and does this occur 
well before the system is in the strict "preformed pair" 
limit? In a related fashion, we address the temperature 
dependence of this breakdown, for a given moderately 
large g value. In this way, the lowest order scheme will 
allow us to establish a criterion for determining T*. 

A. T Matrix: From Weak to Strong Coupling 

The lowest order T-matrix can be written as 

The imaginary component of the electronic self energy 
may be directly computed via 

Im Ek^;^ = ^ [b{ui + eq_k) + /(fq-k)] 
q 

X/miq^^+,^_^(^2_q^2- (6) 

where f{x),b{x) = {e^^ ± (and the real part ob- 

tained via Kramers Kronig transforms). This expression, 
which is not restricted to the lowest order theory, will be 
used throughout the paper. 

It can be seen from Eq. (||) that the self energy is deter- 
mined by Im tq,o. It is instructive to plot this quantity in 
the lowest order theory, as shown in Figure |^, for a range 
of coupling constants and at fixed relative temperature 
T/Tc — 1.5, with wave-vector q = . Here Tc is obtained 
from the usual Thouless criterion for the divergence in 
the T matrix i'^Iq uj=o- The inset shows the full analytic 
T matrix function on the real axis over a more extended 
frequency range and for fixed coupling g/gc = 1-2 and 
T/Tc = 1.1. It can be seen from the main figure that 
the major effect of increased coupling is to introduce a fi- 
nite frequency peak structure into Im iq q whose strength 
grows with g/gc- Using Eq. (^), it follows that this peak 
structure must reflect itself in the self energy. However, 
on the basis of these plots there is no obvious dividing 
line associated with a qualitative change in ImS. Such 



a qualitative change, which signals the breakdown of the 
Fermi liquid state, is expected, somewhere between the 
regime of weak and strong coupling. 

To identify more precisely what is the nature of those 
effects in the T-matrix which are responsible for the 
Fermi liquid breakdown, it is instructive to plot the 
real and imaginary parts of the inverse function tq^j ^ 
which is shown in Figure ^ for moderately strong cou- 
pling g/gc = 1-2, and temperature T/Tc — 1-1. Precisely 
at Tc both the real and imaginary parts of the inverse T 
matrix at q = 0,a; = are zero. However at T/Tc = 1-1 
and for q = it can be seen that the real part of the in- 
verse T matrix is zero at finite frequency, corresponding 
to a pair resonance. Moreover, this resonance disperses 
with wave vector q. It should be noted that there is 
a second zero crossing of the real part of the inverse T- 
matrix which is not associated with a resonance, since for 
this higher frequency the imaginary part of the inverse 
T-matrix is large. 

As shown in the upper right inset, for temperatures 
T/Tc ranging from 1.0 to 1.1, (and for the case g/gc — 
1.0), as the temperature is raised, the resonance gradu- 
ally disappears. This behavior is qualitatively different 
for weaker coupling ( g/gc = 0.6) as shown in the lower 
left inset. In this latter case the only time the real part 
of the inverse T-matrix is zero is in the critical regime, 
which essentially coincides with Tc . The characteristic di- 
viding line between resonant and non-resonant behavior 
appears to be around g/gc ~ 1-0, for this value of the in- 
teraction range. For slightly different parameterizations 
of the interaction this number may change slightly. We 
will show in the next subsection that these resonance 
effects lead to non Fermi liquid like behavior in the elec- 
tronic self energy. 

B. Self Energy: From Weak to Strong Coupling 

To confirm our numerical results, we will compute the 
self energy S^^^ in two different ways. Im ^i^.uj is de- 
rived by direct evaluation of the real frequency integral 
and the Kramers-Kronig relation used to obtain Re Sk,i^. 
We refer to this as the direct real-axis calculation. Al- 
ternatively, we consider the self-energy at the Matsubara 
points using 

q J-oo TT 52 - lUJl - Ek-q 

It should be noted that this expression is equivalent to 
ignoring the Fermi term /(cq-k) in Eq- @- This term is 
not important in the well established pseudogap phase, 
but it is necessary to recover Fermi liquid behavior. A 
Fade approximant scheme |Q is then used to analyti- 
cally continue the self-energy to the real frequency axis. 
These two approaches - the direct real-axis calculation, 
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and the Matsubara scheme - gave rather similar results 
away from the Fermi liquid regime, although only the 
former led to a proper Fermi liquid signature in the self 
energy. The second of these two schemes was used in 
earlier work by our group |22[ . 

In Figure |^ we plot the complex self energy func- 
tion obtained from the Fade approximant approach for 
g/dc = 1-2. For these intermediate coupling strengths, 
a resonance condition in the T-matrix leads to a peak 
in ~Im Sk.w at w -|- eic ~ 0. This is in contrast to the 
Fermi liquid regime where —Im Sk.uj exhibits a minimum 
rather than a maximum. This peak reflects that in the T- 
matrix, however, further amplifled by the zero frequency 
divergent contribution of the Bose factor. These observa- 
tions are confirmed by the real-axis calculations as well; 
the two results for ~Im Sk,w are compared in the inset 
and agree reasonably well. 

The evident association between a peak in —ImY!^^^ 
and a resonance condition in the corresponding T matrix 
(^q^a ^ ) suggests that it is primarily the resonance struc- 
ture ( at low q and low ) in the latter which determines 
the behavior of the self energy. To check this conjecture, 
and to provide some simplification for future calculations, 
in Figure |5|we compare — ImS^^^ deduced using various 
Taylor series expansions in q and U, for the inverse T- 
matrix . The result obtained from the full expression 
is indicated for comparison by the solid curve [Q. 
Reasonable agreement is obtained provided a sufficient 
number of terms (in 17) are included in this Taylor series 
or time dependent Ginzburg- Landau (TDGL) expansion 
of the T- matrix 

<?r' = ^0 - (ao + + bo^^ + C£g9' • (8) 

Here the parameters ag, Qq, bo, and ^^g calculated 
by expanding the full expression for i^^^j . Because of this 
reasonably good agreement, and because of the numerical 
simplicity it introduces, in the remainder of this section 
the direct real axis self energy calculations are based on 
this Landau- Ginzburg form. For our Fade calculations, 
we included the full fi, q dependences. 

It should be stressed that the calculations in the previ- 
ous subsection demonstrate that the full T-matrix t'^]^ 
has a complicated analytic structure. Thus it should 
never be assumed that this function is equivalent to t^'fi- 

Nevertheless, the above calculations of — indicate 
that, for the purposes of computing the self energy, it is 
adequate to introduce a simplified form for the T-matrix. 

The quadratic in fl terms of Eq. (||) were found to 
be reasonably important for obtaining agreement in the 
tails of the self energy function. These quadratic terms 
capture the physics of the second zero crossing of the 
inverse T-matrix, shown in the upper right hand inset of 
Figure |^. We will see below that very near Tc most of 
the physics is dominated by the low frequency regime so 
that iq Q is well approximated by a simple pole structure 
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where Oq = (ro -I- ^l(.c^) / a'^^, and Fq = VLc^{a'^J a'^) . 
Eq. (H), thus, reflects that component of the T-matrix 
which dominates the self energy and related quantities 
near Tc- 

We end this subsection by studying the evolution with 
temperature of — /mE|^°^ as a function of lo for g/gc — 
1.2. The results are plotted in Figure ^ where the three 
curves are at different temperatures (given respectively 
by 1.05, 1.1 and 1.4 T jT^) and correspond to tempera- 
tures just below, at and significantly above the onset for 
pair resonance, which occurs at approximately 1.2rc- A 
Fermi liquid like minimum is present only for the high 
temperature regime. This behavior mirrors the evolu- 
tion with coupling constant seen at fixed T. It should be 
stressed that, upon increasing T, when a single peaked 
spectral function is first recovered, this peak is still quite 
broad. This phase should, thus, not be associated with 
a canonical Fermi liquid. 

C. Parameterization of E and Spectral Function Gap 

In this subsection we identify a pseudogap parame- 
ter and associate this parameter with the calculated gap 
in the spectral function. In the process we arrive at a 
convenient parameterization of the self energy which will 
be useful later on. This parameterization works better 
as temperature decreases towards the well established 
pseudogap regime. Thus this parameterization forms the 
foundation for later work within the full mode coupling 
scheme. 

In Figure ^a) we re-plot the self energy of Figure ^ 
along with the straight line curve for lo. The points at 
which the line intersects —Re Sk,^ are important because 
they lead to structure in the Green's function Gk.^. In 
this way we determine where the real part of the inverse 
Green's function, Re-G^^ = w — ek — -ReSk,^; = 0, cor- 
responding to a set of three zeros. Since the imaginary 
part is small for the two outer roots, the spectral func- 
tion ^k.w, which is given by ^k(^) = — 7r~"'^Im \G\SyiL^ii — > 
LO \- iO)], acquires two peaks separated by 2|Apg|(/5k with 

Note that within our sign conventions |Apgp is a posi- 
tive definite quantity. The spectral function is plotted 
in Figure |^(b). The gapped structure thus corresponds 
to a pseudogap in the electronic properties ( such as 
the density of states). This should be contrasted with 
the single peak in A\^^^ associated with the Fermi liq- 
uid phase (at weaker coupling or higher temperatures). 
The asymmetric broadening of the two peaks is a generic 
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feature of this pseudogap model and is due to the in- 
teraction of correlated pairs with the occupied states 
of the Fermi sea. The dispersion of the two peaks has 
been plotted elsewhere As the momentum vector 

k passes through the Fermi surface the spectral weight 
shifts from the negative to the positive frequency peak. 
Close to the Fermi momentum the peaks disperse roughly 

as = ±^4+ lApfflVk- 

Given the characteristic peak in Im Sk.w , it is useful to 
consider a two parameter model self energy of the form 



^model 



(11) 



A fit of ImT,]^^^ to Eq. ( |ll| ) is presented in Figure g|(a). 
Here the adjusted parameter is 7 since Ap^ can be eval- 
uated from Eq. (^0|). The major difference between the 
two curves is evident primarily in the tails of the func- 
tions. As the temperature gets closer to Tc, we will see 
that the model becomes progressively better. 

In order to further assess the utility of the model, it 
is instructive to examine the dispersion of the peak of 
/mEk^cj. Figure ||(b) plots /m Ek,w for a series of wave 
vectors. This peak dispersion can be compared with that 
of the model, as is shown in the inset of Figure ||(a). 
Agreement is reasonable, and on this basis we will make 
use of Eq. (RtI) in later studies. 



IV. MODE COUPLING THEORY 

In this section we address the three coupled integral 
equations shown in Eqs. (P in the full mode coupled 
regime. These equations must be solved numerically. We 
have just seen that the numerics are tractable for temper- 
atures at and above T*, where the "lowest order" theory 
is adequate. The transition from the Fermi liquid break- 
down to the well established pseudogap phase is difficult 
to treat and its detailed characterization will be deferred 
to a future paper. From the well-established pseudogap 
state, down to T^, we will show that the numerics again 
become tractable. We first address the T-matrix and 
self energy in this regime, concentrating on moderately 
strong coupling {g/gc = 1-2), where there is an apprecia- 
ble pseudogap. We next use these results to determine 
the behavior of Tc as a function of arbitrary g. This lat- 
ter quantity can be derived from the Thouless criterion or 
T-matrix divergence. Mode coupling effects enter in an 
important way, since the T-matrix differs from that de- 
duced in the lowest order theory. The superconductivity 
that results within the mode coupling scheme should be 
viewed as deriving from pairing of pseudogapped elec- 
trons. This superconductivity is, thus, not associated 
with a Fermi liquid state. 



A. Behavior of E and the T-matrix near Tc 

The coupled equations for Sk.w and i"^, as a func- 
tion of temperature, were solved numerically using the 
model self-energy of Eq. ( |Tl| ) . In a given iteration Apg 
and 7 were used to compute iq^n via Eq. (lb). The re- 
sult was then substituted into Eq. (la). The output of 



the latter was then used to extract Apg and 7 via (least 



squares) fits to Eq. (|ll|) and the procedure repeated until 
convergence. Once self-consistency in these parameters 
is reached, the resulting self-energy and T-matrix can 
be calculated as an output from the full equations ( p^ 
lb). To distinguish these from the modeled functional 
forms, in the following we refer to these results as the 
full functions for Imtq^Q and /mSk,i^. Unless indicated 
otherwise, these are the quantities plotted throughout 
the remainder of this paper. 

The resulting inverse T- matrix is plotted in Fig- 
ure ^ The upper right inset plots Im o and the 
lower inset — /mEk,^ for a temperature slightly above 
Tc- (Proximity to Tc is signaled by an incipient diver- 
gence in the T-matrix). Here, g/gc = 1.2. This behavior 
should be compared with the lowest order results in Fig- 
ures HjIIJ. The salient feature in the main body of the 
figure is that a resonance appears, just as in the lowest 
order calculations. What is striking, however, is that the 
resonance is sharper as a consequence of feedback effects, 
due to a suppressed, essentially gapped, imaginary part. 
The (approximate) gap in the imaginary part is to be 
directly associated with the pseudogap in the electronic 
spectrum. Thus we have the interesting result that feed- 
back effects stabilize the resonance found in the lowest 
order theory. Stated in more physical terms, as Tc is ap- 
proached, pairs of low momenta become very long lived. 

The upper right inset shows Im tq_n which exhibits a 
significantly sharper peak structure than was seen in the 
lowest order theory. It is to be expected that this peak 
will reflect itself in — /mEk.w, which is plotted in the lower 
left inset. It should be noted that the figure plots the 
full function not that of Eq. ( pi] ) to which it can be 
well fitted. Indeed the fit becomes progressively better 
as temperature is lowered, and concomitantly 7 becomes 
smaller. 

This analysis will be used later to determine the ac- 
tual temperature at which the T-matrix diverges. As the 
temperature is lowered towards Tc several important sim- 
plifications occur: (i) 7 approaches zero and (ii) most of 
the contribution to the self energy is captured by approx- 
imating the T-matrix inside the integral of Eq. (^ by the 
TDGL form. Point (i) indicates that the self energy be- 
havior near Tc is a continuation of the (much weaker) 
peak in Im Sk.i.^ at w = — ek which was found earlier as 
the Fermi liquid breaks down near T*. With lower T 
this peak grows progressively sharper until the instabil- 
ity at Tc is reached. The TDGL expansion noted in (ii) 
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is given by Eq. (g). This equation reflects the known 
analytic form (at low 17 and q) of the divergence in the 
T-matrix. It is not unexpected that this incipient diver- 
gence will dominate the properties of the electronic self 
energy. Moreover, this simplification is readily incorpo- 
rated into the numerical iteration scheme discussed above 
to expedite our calculations. As an intermediate step, we 
fit the calculated T-matrix to Eq. (0) for use in Eq. (la). 



It should be stressed that after convergence the full self 
energy and T- matrix are produced. Moreover, at each 
step the validity of Eq. ( pi] ) and Eq. (||) are checked 
by comparing with their respective counterparts given in 
Eq. (0) and Eq. (0). 

We conclude this subsection by comparing the evolu- 
tion of the single particle parameters Apg and 7 with 
their pair counterparts which we parameterize as Qq and 
Fq by fitting the resonance feature in the T-matrix to 
Eq. (^. The pair parameters, ilq and Fq, thus, corre- 
spond respectively to the energy and width of the res- 
onance in the T-matrix. The results of our numerical 
iteration scheme are shown in Figure |lO| f or g/gc = 1-2. 
The shaded regions indicate where Eq. (|ll[) breaks down, 
as does the lowest order theory. Even though the pseudo- 
gap persists, the system begins to cross-over towards the 
Fermi liquid. Thus the parameters Apg and 7 are indi- 
cated only for a limited range of temperatures, where Eq. 
( pi] ) is relevant. The temperature for Fermi liquid onset 
~ T*, is just beyond the shaded region corresponding 
to the regime of validity of the lowest order theory. At 
higher temperatures Fq^o and r2q=o are plotted to join, 
after extrapolation, onto the results of the lowest order 
theory. The temperature T* can be read off from the 
endpoint of the two curves in Figure ^(b). This is the 
highest temperature at which F and are defined. 

This figure re-enforces the earlier observation that the 
pairs are relatively long lived over much of this temper- 
ature regime: the decay rate, Fq^o, remains comparable 
to riq=o- The small value of Fq=o is also reflected in 
the small value of 7. We may view this (well established 
pseudogap phase, to the left of the shaded region in Fig- 
ure p^ as an extended critical regime which arises as a 
result of the gapping of the T-matrix via feedback effects. 

This variation in 7, and Apg, along with our earlier 
analysis from the lowest order theory, (see, for example, 
Figure ^) has direct implications for experiment: the 
spectral function A^j^ exhibits a single broad peak at 
T > T*, which evolves into two peaks at T « T*. The 
separation between these peaks increases with decreas- 
ing temperature, as Apg increases. At the same time 
since 7 decreases as Tc is approached, the two peaks pro- 
gressively sharpen in the vicinity of the superconduct- 
ing transition. This behavior represents, then, a smooth 
transition of pseudogap behavior into the superconduct- 
ing state. It should not, however, be assumed that the 
pseudogap parameter Apg is equivalent to the supercon- 
ducting gap Age below Tc. This latter issue will be dis- 



cussed in a future paper 134]. 



B. Calculation of Tc 

In this subsection we deduce the transition temper- 
ature for variable coupling constant using two different 
numerical approaches. We use the iterative numerical 
scheme discussed above to determine Tc by studying the 
diverging T-matrix as the temperature is lowered. In 
this way Tc is obtained as the temperature at which both 
Fq^o, and ilq=o are identically zero. 

This scheme yields numerically equivalent answers to 
those obtained using an alternative and simpler ap- 
proach. For the purposes of calculating Tc, the largest 
contribution to the self energy in Eq. (|la|), may be seen 
to come from the low frequency, long wavelength phase 
space region where the T-matrix is large. The integral is 
well approximated by |}45| 



A 2 2 ^(0) 



(12) 



Indeed, this equation is equivalent to Eq. (11) in which 
7 is set to zero, as was found to be the case numerically 
at Tc. 

It follows from this equation that 



■-0,0 



-2/(gk) „2 



(13) 



where i?k = y^k + ^pg'^'k- pseudogap parameter 

which appears in Eq. ( p^ ) is given by 



Ai 



E 



(14) 



The lowest order analogue of this equation was inferred 
earlier in Section 3. It is evident that Apg in Eq. ( p^ ) 
coincides with the square amplitude of pairing fluctua- 
tions, g'^ {c' c'^ cc) , which can be defined more generally 
away from Tc. 

In order to evaluate Tc, Eq. ( p^ must be combined 
with the Thouless criterion from Eq. (|l3|), along with 
the number equation. The three fundamental equations 
for Tc are Eq. dlj) , along with 



1 - 2/(£;k) 2 n 



2Eu 



2E 



,,2 I iii 



(15) 
(16) 



where = (1 



It should be noted that Eqs. (|l|), (|l|) and (^ show 
that the single particle Green's function, the Thouless 
criterion and the number equation assume a particularly 
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simple form which coincides which their (below Tc) coun- 
terparts obtained in the standard superconducting theory 
but with the pseudogap playing the role of the super- 
conducting gap. This simplicity (as well as the related 
BCS-like behavior, at small g, below Tc) would not 
obtain if fully rcnormalized Green's functions are used 
everywhere ||, [0- Eqs.(|l|)-(|l|) can be viewed 
as a simple (one variable, Apg) parameterization of Eqs. 
(0), (0), and dl^. 

The overall behavior of Tc is compared with that ob- 
tained from the NSR approximation of Ref . |^ , as well 
as that of strict BCS theory in Figure |ll|. The non- 
monotonic behavior of Tc results from the facts that (1) 
the high g asymptote must be approached from below 
1^ and (2) as in Ref. |2^, the low g exponential de- 
pendence tends to overshoot this asymptote. Here, the 
overshoot is even more marked than in the NSR theory 
of Ref because of self energy effects which pin Ep at 
/i, as can be seen from the inset. Thus, at small g our 
mode coupling curves tend to follow the BCS result more 
closely. Point (1) is an important point. As g increases, 
the pair size decreases and consequently the Pauli princi- 
ple repulsion is diminished. In this way Tc must increase 
with g at large g in contrast to the behavior found in Ref 

ii- 

Figure ^ presents a comparison of the three different 
energy scales: the chemical potential fi, the pseudogap 
parameter Apg and Tc- It can be seen that the maxi- 
mum in Tc is associated with fi ~ Apg and the mini- 
mum with /i = 0. Indeed, the complex behavior of Tc 
shown in Figure |l^ can be understood on general phys- 
ical grounds. A local maximum appears in the Tc curve 
as a consequence of a growing (with increased coupling) 
pseudogap Apg in the fermionic spectrum which weak- 
ens the superconductivity. However, even as Apg grows, 
superconductivity is sustained. In the present scenario, 
superconductivity is preserved by the conversion of an 
increasing fraction of fermions to bosonic states, which 
can then Bose condense. Once the fermionic conversion 
is complete (/i — 0), Tc begins to increase again with 
coupling. 

The behavior of Tc on an expanded coupling constant 
scale, for different ranges of the interaction (parameter- 
ized by ko/kp, which enters via (fk = {1 + k^ /kQ)~^^'^) 
is shown in Figure The limiting value of Tc for large 
values of g/gc approaches the ideal Bose-Einstein con- 
densation temperature Tbe = 0.218£'f for the case of 
short range attraction so that /cq ^ oo. The qualita- 
tive shape of the Tc curve, however, is retained for ko/k-p 
greater than about 0.5. For larger range interactions, the 
solution disappears for some g/gc', then, when /i is suffi- 
ciently negative, the transition reappears approaching a 
continuously (with kg/k-p 0) decreasing asymptote. 



C. Phase Diagram 

Our results for Tc as a function oi g/gc can be consoli- 
dated into a quantitative phase diagram shown in Figure 
which includes only the physical parameter regime 
(/i > 0). We, thus, do not show larger values of the 
coupling constant where the system acquires a bosonic 
or pre-formed pair character. The parameter T* is also 
indicated for each coupling constant. This was deter- 
mined as the temperature at which a resonance at q = 
first appears in the T-matrix. For lower values of the 
coupling constant g/gc < 1-0, this resonance onset essen- 
tially coincides with Tc, but above this regime there is a 
clear separation between the two temperature scales. For 
g/dc ~ 1-5 , as the pseudogap size further increases, it 
starts to weaken the superconductivity, in large part be- 
cause of the decreased density of fermionic states which 
can participate in the pairing. 

There is, thus, a delicate balance between supercon- 
ducting pairing and the opening of a pseudogap associ- 
ated with strong superconducting attraction. In general 
gaps in the spectral function, as in charge or spin density 
wave systems are expected to weaken superconductiv- 
ity. What is distinctive about a gap which arises from 
the pairing, is that superconductivity persists even when 
the gap is an order of magnitude greater than Tc- This 
occurs in the present model (as seen in Figure [ij) and 
this also occurs experimentally in the cuprates. Su- 
perconductivity can coexist with large spectral function 
gaps when these gaps are associated with a transition to 
a bosonic state: once /j, becomes negative with increasing 
g, the conversion of fermions to bosons is complete and 
superconductivity, through Bose- Einstein condensation, 
occurs in an unhindered fashion. We view the experimen- 
tally observed large size of Apg /Tc as suggestive of 
pseudogaps associated with superconducting pairing. It 
is doubtful that other pseudogap mechanisms would be 
as compatible with superconductivity. 

Our phase diagram is plotted as a function g/gc- The 
usual cuprate phase diagram is plotted in terms of hole 
concentration dependence x. This latter quantity has not 
been discussed thus far. Nevertheless, there are two im- 
portant effects which are associated with the transition 
to the insulator a: ^ 0, which both enhance pseudogap 
effects with decreasing x. Here we assume for definite- 
ness that the pairing mechanism is not changed as the 
insulator is approached. (This paper has not made any 
presumptions about this mechanism and any associated 
hole concentration dependence is not relevant to our con- 
siderations) . As the insulator is approached, the system 
generally becomes progressively more two dimensional, 
as seen for example, by the c-axis resistivity |46|. The 
transition to two dimensionality is associated [ [40| with 
smaller and smaller Tc, so that the size of the pseudogap 
temperature regime, which varies as T* /Tc, will grow as 
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the insulator is approached. 

Another important component of the transition to the 
insulator is the decrease in plasma frequency. This is- 
sue has been emphasized by Emery and Kivelson p2[ . 
Earlier work by our group | |25[ | demonstrated that when 
the Coulomb channel, along with charge effects, were in- 
cluded into the Nozieres Schmitt-Rink formulation, they 
helped to stabilize the mean field regime. The smaller 
the plasma frequency the greater the deviation from BCS 
behavior. This observation is not readily associated with 
softening of phase fluctuations in the sense of Ref. . 
In the present approach, the Coulomb channel creates 
a particle-hole attraction which is competitive with the 
particle particle (Cooper pair) attraction. In this way, 
one may infer that as the insulator is approached, the 
system should become progressively less mean field like, 
so that pseudogap effects are expected to be enhanced. 
These and related effects will be discussed in more detail 
in a companion paper po|. 



D. Relation to the literature 

In the previous section we discussed self energy effects, 
Tc calculations and, in passing, some aspects of time de- 
pendent Ginzburg-Landau (TDGL) theory which should 
be compared to related calculations in the literature. 

Our self energy calculations are in many ways similar 
to related work on the effects on electronic properties of 
proximity to magnetic transitions [ ^7p^ . By contrast, 
in these other systems, there have been no reports of a 
breakdown of the Fermi liquid phase, except in the crit- 
ical regime. In most of these previous calculations, as 
in some of ours, the relevant susceptibility was expanded 
around the long wave length, low frequency limit. What 
causes the Fermi liquid breakdown in our case, can be 
traced to a special feature in this form of the pair sus- 
ceptibility, which has no counterpart in these other prob- 
lems. As noted earlier, such a breakdown is naturally ex- 
pected in any interpolation scheme which varies smoothly 
between the fermionic and bosonic limits of g/gc- Our 
calculations show that the Fermi liquid like character in 
the self energy is found to disappear as resonant struc- 
ture in the T-matrix sets in. This can be seen directly 
from evaluating the integral Eq. (||) using Eq. (||). Here 
resonance effects enter via the ratio (oQ/aQ). The larger 
is this parameter, the more pronounced are pair reso- 
nances. Indeed, we are able to tune from Fermi liquid to 
non Fermi liquid behavior simply by increasing the size 
of this ratio, in very much the same way as was seen in 
the full calculations, by increasing g/gc- For values of 
the ratio near « 1 the pseudogap state appears. Mode 
coupling effects which renormalize (aQ/ag), can also be 
incorporated at this TDGL level of approximation, as 
noted below. 



Our TDGL expansion of the T-matrix, Eq. (g) can be 
compared to that discussed in the context of earlier work 
on the BCS Bose-Einstein cross-over problem It 
is well known that as the Bose Einstein regime is ap- 
proached, the parameter Oq increases relative to ag. In 
this way the T-matrix acquires a propagating component 
in addition to the diffusive term of BCS theory. Indeed, 
this is precisely the situation in which the T matrix is 
described by a resonance. What is new in the present 
theory is that, not only do we explicitly associate this res- 
onance with the break down of the Fermi liquid, but also 
for our mode coupling calculations, the diffusive compo- 
nent (e.g., the parameter Oq) remains small for an ex- 
tended range of temperatures above Tc and essentially 
vanishes at Tc for all values oi g/gc- This represents a re- 
statement of our earlier conclusions that mode coupling 
effects lead to extremely long lived pairs. These long lived 
pairs are a reflection of the long lived single particle states 
whose rate of decay is characterized by 7. Indeed, it can 
be shown by direct expansion of the calculated T-matrix 
that 



7V(0). 



A |2 ■ 



(17) 



so that Oq varies as a power of (T — Tc) sufficiently close 
to Tc- This should be contrasted with the result obtained 
p^ , ^ without mode coupling effects. 

Finally, in the two dimensional (2d) limit our Tc cal- 
culations pO| can be compared with related work in the 
literature. This 2d limit was noted to be problematic fol- 
lowing earlier work by Schmitt-Rink, Varma and Ruck- 
enstein [Q, who applied the NSR theory to the low di- 
mensional case. It was found that, while Tc was zero 
as expected, ^ was also negative for all values of g/gc, 
including also in the weak coupling limit. Serene |371 ] 
pointed out this latter unphysical result arose from the 
use of the truncated Dyson equation of Eq. (||) . Yamada 
and co-workers ||5l|] attempted to correct this problem by 
including the lowest order "box" diagram in the T-matrix 
(See Figure ^(c)). In this way some mode coupling effects 
were included; it was found that /i « Ep at arbitrarily 
weak coupling, while was zero, as expected in the 2d 
limit. However, the authors were unable to find a the- 
ory which smoothly interpolated between the weak and 
strong coupling limits. 

The present calculations (which, in contrast to earlier 
"mode-coupling schemes" of Ref. |^,|3|,|l] , include a se- 
ries of higher order diagrams ||5^), yield Tc = 0, with /i 
varying continuously from Ep to the large negative val- 
ues of the strong coupling limit. The vanishing of Tc 
arises via Eq. (14), which because of 2d phase space re- 
strictions on the Bose factor, yields an arbitrarily large 
Apg which makes it impossible to satisfy the Thouless 
criterion. In this way, problems with earlier 2d crossover 
theories are corrected pO[ . 

We end this discussion by noting that non monotonic 
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behavior in Tc with coupling g has ample precedent in 
more conventional Fermi liquid-based schemes, where it 
is found that, in some cases, may decrease with in- 
creased large coupling constant, or saturate with g. This 
follows within Eliashberg theory as a consequence of the 
competing mass enhancement factor. The resulting form 
for Tc arises from a competition between the pairing cou- 
pling constant and the mass enhancement ratio ]55| ] . The 
present calculations report a similar effect, which arises, 
however, in the non Fermi liquid regime. Here the com- 
petition is represented by the pseudogap parameter Apg , 
which like m* /m, can undermine the effectiveness of the 
superconducting attraction. 

V. CONCLUSIONS 

In this paper we have applied a diagrammatic decou- 
pling scheme proposed by Kadanoff and Martin, (which 
was further extended by Patton), to the BCS- Bosc Ein- 
stein cross over problem. The strongest support for this 
theoretical framework is that it reproduces the BCS limit 
at very weak coupling, in contrast to alternate diagram- 
matic approaches in the literature. The results of our 
extensive numerical calculations are summarized in Fig- 
ure in which we plot our schematic physical picture, 
as well as the imaginary component of the T-matrix, the 
imaginary part of the self energy and spectral function 
for each of the three regimes: the weak coupling, BCS 
or "Fermi liquid" case, the "pseudogap state" and the 
"bosonic state", or pre- formed pair limit. The physical 
picture of each respective normal state corresponds to 
free fermions, resonantly scattered fermions and bound 
fermions. The curves shown in the figure for the first two 
cases represent quantitative plots, related to or equiva- 
lent to those we have presented throughout this paper. 

The present calculations have addressed a three dimen- 
sional s- wave jellium model and have not made direct 
contact with the cuprate superconductors. The effects of 
introducing d-wave pairing on a quasi two dimensional 
lattice will be discussed elsewhere these effects do 
not qualitatively undermine our conclusions. Indeed, at 
this stage we can make preliminary contact with exper- 
iments. In the pseudogap phase, our spectral functions 
are generally consistent with a leading edge (albeit s- 
wave) gap in the photoeniission spectra. The gap, Apg, 
grows in magnitude with decreasing temperature, while 
at the same time there is a narrowing of the associated 
spectral function peak, which enters via the parameter 7. 
Precisely at Tc, 7 is found to vanish, as appears to be the 
case experimentally [S^ ]. Away from Tc it grows rather 
slowly with increasing temperature, presumably because 
of the assumption of s- wave pairing. 

Our calculations have revealed a rich structure in the 
phase diagram (See Figure ^4|) for T* , and Tc as a func- 
tion of coupling constant. It is most likely that the 



cuprates are in the regime of moderate coupling so that 
the pseudogap is not particularly large compared to elec- 
tronic energies [^[|^. Combined with other experimental 
constraints, this appears consistent with our model only 
in the quasi 2d limit |Q. However, it must be stressed 
that, because it makes no assumption about the pair- 
ing mechanism, the present theory is not appropriate for 
addressing the detailed hole concentration dependence of 
the cuprate phase diagram, except in this very qualitative 
fashion. A more detailed discussion of this and related 
issues is presented in a companion paper . 

Our diagrammatic calculations of the cross over from 
BCS to Bose Einstein behavior, with increased coupling 
have several satisfying features . Among these, we find 
(1) a well established pseudogap phase with a spectral 
function which evolves continuously into the characteris- 
tic two peaked structure of a conventional superconduc- 
tor, (2) that the Bose Einstein asymptote of our Tc curves 
is properly approached from below, (3) that a BCS limit 
is embedded in our weak coupling expression for Tc, and 
(4) that, as will be discussed in more detail in a subse- 
quent paper, the 2d (and quasi 2d ) limits of our theory 
appear to be well behaved, smoothly varying with cou- 
pling, and yield a reasonable parameterization of cuprate 
energy scales. Consequently, at the present time, we be- 
lieve that these cross-over theories, in the regime of inter- 
mediate coupling, should be viewed as viable candidates 
for characterizing the precursor superconductivity con- 
tributions to the cuprate pseudogap. 
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FIG. 1. The class of diagrams considered in this paper for 
the T-matrix and self-energy: (a) The pairing approxima- 
tion; (b) the lowest order theory; (c) Dyson's equation for the 
T-matrix. 




FIG. 2. Main figure: Frequency dependence of Imt^J^Q q 
for T/Tc — 1.5 and coupling (a) g/gc ~ 1-5; {h)g/gc — 1.2; 
(c) g/gc = 1.0; (d) g/gc = 0.6. Inset: Real (solid) and imag- 
inary (dashed) components of T matrix over extended fre- 
quency range, with g/gc ~ 1.2 and T/Tc = 1.1. In this and all 
subsequent figures energies are measured in units of Ep. 




FIG. 3. Main figure: Resonance efi'ects illustrated by 
real (solid) and imaginary parts (dashed line) of piq'HoLj for 
g/gc = 1.2, and T/Tc = 1.1. Upper inset: Regt'-°1~^^ (solid 
curves) for T/Tc from 1.0 (bottom curve) to 1.1 (top curve), 
and g/gc = 1.2. Dashed curve is Imgt^lg]^, which is un- 
changed for these T/Tc. Lower inset: Same as in upper 
inset with g/gc = 0.6. 
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FIG. 4. Main figure: ReEkp,[j (solid line), and ImEkp,[j 
(dashed line), within the lowest order theory. Inset shows a 
comparison between results for ImEkp.ij using direct real axis 
(solid ) and Matsubara (dashed line) schemes. 
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FIG. 5. Comparison of the self-energy (heavy line), 
with various level approximations (see legend) based on ex- 



panding gtq}j 



at low q, Q. 
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FIG. 6. Temperature dependence of 
tion of tj, for T/Tc = 1.05, 1.1 and 1.4. As the temperature 
is lowered, the quadratic minimum at zero frequency - the 
sahent feature of a Landau Fermi liquid - is gradually turned 
into a maximum in the pseudogap regime. 
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FIG. 7. (a) Frequency dependence of E^'^ ^. The single 
particle spectral function will have a peak (at u indicated 
by the circles) whenever ReEkp.t^ crosses the unit diagonal 
(shown by the line) and ImEkp.cj is small, (b) Associated 
spectral function j4k,t^ obtained from the self-energy plotted 
in (a). 
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FIG. 8. (a) Main figure: Numerically calculated 
— /rra Ek,ij (heavy line) and Lorentzian fit (dashed line) pa- 
rameterized by Eq. (|ll|). Inset: The dispersion of peak 
in full (solid line) and model (dashed line) self energy; (b) 
— /m Ek.ij for k/kp = 0.1 (bottom curve) to k/kp ~ 1.5 (top 
curve) . The curves are vertically offset for clarity. 
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FIG. 9. Main figure: The real (solid) and imaginary part 
(dashed line) of inverse T-matrix gt~\^ at full mode coupling 
level, for temperatures slightly above Tc- Upper riglit inset 
plots Imtci,n and Lower left inset plots — /mEk.i^ at the 
same temperature. 
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FIG. 10. The evolution of the single particle parameters 
Apg and 7 (a), contrasted with the pair parameters f2q and 
Fq (b). The coupling is fixed at g/gc = 1-2. For temperatures 
above and below the shaded area, the lowest order theory and 
mode coupling scheme are used, respectively. Curves inside 
the shaded region were obtained by interpolation. 




FIG. 11. Comparison of the Tc curves within the BCS the- 
ory (dotted line), the Nozieres-Schmitt-Rink scheme (dashed 
line), and our mode coupling results (heavy line). The inset 

shows the dependence of the chemical potential n{T = Tc,g) 
on the coupling constant for these three cases. 
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FIG. 12. The evolution of Tc, n, and |A| with coupling 
strength in self-consistent mode coupling theory. 
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FIG. 13. Tc versus coupling constant over an extended 
range for three different values of the parameter kp/ko- The 
minimum in the two larger ko curves is associated with = 0. 
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FIG. 14. Quantitative phase diagram in the fermionic 

regime, /i > 0. The shaded region indicates whore there is 
a well established pseudogap in the single particle spectral 
function. 
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FIG. 15. Summary of our conclusions: schematic plot of 
the resonant pair or pseudogap regime as a crossover between 
free fermions (BCS) and bound fermions (Bose-Einstein) with 
fermions indicated by the solid dots, dashed lines around pairs 
denote mctastablc pair states, while stable preformed pairs 
arc represented by encircled pairs. Panels (a), (b) and (c) 
illustrate the evolution with coupling of: (a) /mtq=o.s), (b) 
— /m.Ek,aj, and (c) ^k,uj- The vertical scales of figures for BCS 
and resonance regime have appeared elsewhere, throughout 
the paper. 
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